Spin polarization in open quantum dots 
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We investigate coherent transport through open lateral quantum dots using recursive Green's 
function technique, incorporating exchange-correlation effects within the Density Functional Theory 
(DFT) in the local spin-density approximation (LSDA) . At low electron densities the current is spin- 
polarized and electron density in the dot shows a strong spin polarization. As the electron density 
increases the spin polarization in the dot gradually diminishes. These findings are consistent with 
available experimental observations. Results of our DFT-based modeling indicate that utilization of 
the simplified approaches that use phenomenological parameters and/or model Hamiltonians might 
not be always reliable for theoretical predictions as well as interpretations of the experiments. 
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Introduction. A detailed understanding of spin-related 
phenomena in quantum systems is necessary for future 
spintronics applications such as spin filtering devices^], 
spin-FETs Q|, nonvolatile computer memories 0, etc. 
Semiconductor quantum wires, dots and antidots de- 
fined in a two-dimensional electron gas (2DEG) represent 
promising systems for the implementation of quantum 
spintronic devices In this context, a topical question 
is whether the spin degeneracy in these structures can be 
lifted. 

Experimental studies indicate the existence of a spon- 
taneous spin polarization at low densities in the 2DEG0. 
The spontaneous spin polarization was suggested as the 
origin of "0.7-structure" in the conductance of a quantum 
points contact (QPC) 0. Concerning spin polarization 
in open quantum dots, i.e., dots with strong coupling to 
leads as opposed to the Coulomb blockade regime, the ex- 
isting experiments show conflicting findings. A statistical 
analysis of conductance fluctuations indicated a spin 
degeneracy at low magnetic fields. The similar conclu- 
sion follows from the results of Folk et al. ja] who experi- 
mentally demonstrated the operation of an open dot as a 
spin filter. In contrast, low-field magnctoconductance of 
small dots Q showed a pronounced peak splitting that 
was taken as a signature of the spin polarization. In 
contrast to QPCs, where theoretical investigations of the 
spin polarization have been a subject of lively discussions 
during recent years [H E 111 III El, theoretical 
description of spin-polarization effects in open quantum 
dots has received far less attention. The main purpose of 
this paper is to provide such the description. 

Modelling transport through quantum dots can be 
done from conceptually different standpoints. For ex- 
ample, charging and coupling constants might be taken 
as phenomenological parameters of the theory within a 
model Hamiltonian j^, 0|- It is however not always 
evident whether such a simplified description is suffi- 
cient to capture the essential physics, and it is not al- 
ways straightforward to relate quantitatively the above 
parameters to the physical processes they represent in 
the real system. Another common approach is to ap- 
proximate the open system at hand (which is described 



by the extended (scattering) states), by a corresponding 
closed system, and then to calculate a self-consistent po- 
tential of this system. The calculated potential is then 
used to obtain the transmission coefficient on the basis of 
e.g. transfer matrix technique or related methods [rR fl8l | . 
This approach avoids ambiguity related to the choice of a 
model Hamiltonian and phenomenological constants but 
depends on the crucial (and not obvious) assumption that 
there is no qualitative difference between many-electron 
effective potentials for open and closed systems. 

In the present paper we present an approach based on 
full quantum mechanical many-electron magnetotrans- 
port calculations for open system. Using the recursive 
Green's function technique [19j we compute the scatter- 
ing solutions of the two-dimensional Schrodinger equa- 
tion in a magnetic field. Following the parametriza- 
tion for the exchange and correlation energy functionals 
of Tanatar and Cerperly [2(|, the electron-electron in- 
teraction is incorporated within the Density Functional 
Theory (DFT) in the Local Spin Density Approxima- 
tion (LSDA). The choice of DFT+LSDA for the descrip- 
tion of many-electron effects is motivated, on one hand, 
by its efficiency in practical implementation within a 
standard Kohn-Sham formalism |2l|. and, on the other 
hand, by the exellent agreement between the DFT and 
exact diagonalization|22j and variational Monte-Carlo 
calculations [H [3] performed for few-electron quantum 
dots . 

We stress that our approach involves no phenomeno- 
logical parameters (like coupling/charging constants), 
and that the self-consistent potential is computed for the 
scattering states of an open system. As a note, conceptu- 
ally similar approaches have been used to model quantum 
transport in molecular conductance and related molecu- 
lar systems |25j as well as to study quantum point con- 
tacts in zero mag netic field [H 03. 

Model. We consider a two dimensional (2D) lateral 
quantum dot defined by a top gate in a rectangular ge- 
ometry. The gates induce an external electrostatic con- 
finement described by a harmonic potential dim 
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FIG. 1: External confinement potential v ext and the computa- 
tional area used for the dot. The strength of the confinement 
is hu x = 0.2meV, hw y — 0.15meV, the area size is 360x210nm 
and the lead width is 80nm. Shaded area in the computational 
domain indicates schematically the region where the confining 
potential lies below the Fermi level. 



The dot is connected to electron reservoirs by two semi- 
infinite leads, Fig. ^ The total confining potential in- 
cludes also a contribution from a positive background 
due to a uniform layer of donors By changing the 

donor concentration p + , which contributes to the classi- 
cal Hartree potential (Eq.EJ, we can effectively alter the 
electron density in the dot. In this manner, the transport 
characteristics of dots containing different electron densi- 
ties can be investigated. In the present paper we present 
a case study of the magnctotransport in a dot with "low" , 
"intermediate" and "high" electron densities. This cor- 
responds to a donor concentration p + /10 14 m -2 = 2, 8, 20 
(resulting in similar average electron densities, n s ). The 
interaction parameter r s = l/a* B ^/rt 



:n s , where a* B is the 
effective Bohr radius, is 3.75, 2.00 and 1.28 respectively. 

The dot is described by the effective Hamiltonian in a 
perpendicular magnetic field, = Bi, 



Hqd = H + <//(r) + gp B Ba, 



where Hq is 



H = - 



2m* 



d eiBy 
dx h 
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In Eq. © the Landau gauge is chosen, A = (—By, 0,0); 
77i* is the effective mass in GaAs (= 0.067m e ). The last 
term in Eq. accounts for Zeeman energy where fig = 
eh/2m e is the Bohr magneton, a = ±^ describes spin- 
up and spin-down states, f, |, and g factor of GaAs is 
g = -0.44. 

The effective potential, v°f j (r) , follows from the Kohn- 
Sham theory |2ll. and consists of three parts, 



<f f W = v H (r) + < c (r) + v ext {r). 



(4) 



vh(f) is the classical Hartree potential due to the elec- 
tron density n(r) = J2a n<T ( r )^ anc ^ the l a y er °f positive 



vh{f) 



47T£o£r 



dt 



,n(r')-p+ 



(5) 



The exchange and correlation potential v xc (r) = v x [rj 
v c (r) in the local spin density approximation reads |21 



vZ c =—(n°e xc (n,C(r)) 



(G) 



where ^ = " T T"j. is the spin-polarization. For e xc we 
have used the paramctrization of Tanatar and Cerperly 
[20j : the explicit expressions for v x (r) and v c (r) can be 
found in [j^ . 

The leads are modeled as hard wall potentials where 
the magnetic field is restricted to zero. This is fully jus- 
tified in the field interval under consideration (B < 0.5 
T), because the distortion of the wave function in the 
leads due to the effect of magnetic field is negligible. All 
the results presented here correspond to one propagating 
mode in the leads. 

Calculations. In order to compute the electron den- 
sity and the conductance of the dot we use the recursive 
Green's function techniaue[lflj. We discretize Eq. (J2J and 
introduce the tight-binding Hamiltonian (with the lat- 
tice constant of a = 10 nm), where the perpendicular 
magnetic field is included by Peierl's substitution [28j |. 
Introducing the retarded Green's function, 



(E - H a 



(7) 



we compute the density at site r through [2^] 

n°{v) = --( Im[g°(r,r,E)]f(E-E F )dE, (8) 

t J -00 

where / is the Fermi-Dirac distribution and Ef is 
the Fermi energy. Equations ©~© are solved self- 
consistently for the spin up/down densities. Note that a 
direct integration along the real axis in Eq. (JSJ) is rather 
ineffective as its numerical accuracy is not sufficient to 
achieve convergence of the self-consistent electron den- 
sity. Because of this, we transform integration contour 
into the complex plane ^s(E) > 0, where the Green's 
function is much more smoother. This also allows us 
to include bound states in the dot, with energies below 
any propagating mode in the leads. We would like to 
stress that the Hamiltonian H a in Eq. J7J) (and thus 
the corresponding Green's function Q a ) correspond to an 
open system representing the quantum dots and semi- 
infinite leads. Using the converged, self-consistent po- 
tential, the transmission coefficient T a is computed from 
the Green's function between the two leads, T a (E) oc 
E m ,» (V>« \G° (rL R ,r Lh , E)\t/j m ), where 

f Lr l denotes the 

coordinate at the right respectively left lead, and ip n ,m 
the wave functions in the leads[2^. Finally, the conduc- 
tance G = ^2 G (j is found from the Landauer formula, 
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FIG. 2: (color online) (a) The density for the spin-down 
electrons n^(r) in the low density dot calculated within the 
Hartree approximation (note that n^(r) and n T (r) are indis- 
tinguishable on the scale of figure) . (b) Average electron den- 
sity for spin-down and spin-up electrons in the dot, nj , vs. 
magnetic field, (c) Conductance (in units of e 2 /h) through 
the dot vs. magnetic field. Temperature T=2K. 
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FIG. 3: (color online) The density for the spin-up and spin- 
down electrons ra T (r) , ra^ (r) (left and right panels respectively) 
in a (a) low, (b) intermediate and (c) high density dot calcu- 
lated within LSD A. T=2K. 

Results. To outline the role of exchange and correla- 
tion interactions we first study the magnetotransport in 
a quantum dot within the Hartree approximation (i.e., 
when v xc (r) is not included in the effective potential (QJ). 
Figure |2Ia),(b) show the electron densities n^(r), n^(r) 
and the average spin-up/down density, n\/ h*, in the 
low density quantum dot. As expected, the electron den- 
sities and the effective potential are practically the same 
for spin- up/down electrons. Although the Zeeman term 
lifts the degeneracy between the electrons of different 
spins, the difference between spin- up/down electrons is 
very small {\n\/n\\ ~ 1.01). 



Figures|3{a) and^Ja) show the electron densities n^(r), 
n^(r), the average electron densities, n], nj;, and the con- 
ductance for the small quantum dot calculated in the 
framework of LSDA, where the effective potential is given 
by Eq. QJ. The difference with the Hartree approxima- 
tion (Fig. [2J is striking. The exchange interaction causes 
a significant spin polarization with |n^/nj| ~ 3. Note 
that spin-up and spin-down electrons are spatially sep- 
arated and localized in different parts of the dot. It is 
interesting to note that one of the spin channels is com- 
pletely suppressed such the conductance is dominated 
by the spin-down electrons. This findings indicate that 
small quantum dots can be used for injection of the spin- 
polarized current. It is also important to stress that the 
spin polarization persists for B = 0, which clearly points 
out to the exchange (not Zeeman) interaction as the ori- 
gin of the effect. 

With the increase of the electron density the spin po- 
larization in the dots decreases. This trend is illus- 
trated in Fig. 01 (b) where the spin polarization in the 
intermediate density dot is decreased to |n^/nj| ~ 1.2, 
and for the high density dot is practically disappeared 
|n^/n]| ~ 1.01. This behavior is also clearly manifested 
in the distribution of the electron densities, Fig. In 
the intermediate density dot the spin- up/down densities 

(r) , (r) are no longer as distinctly separated as in the 
small one, and in the large dot n^(r), n^(r) are, except at 
the edges, practically the same throughout the dot. The 
peculiarities at the edges in^c) is due to the finite size of 
our computational domain; as the density and Coulomb 
interaction increase the electrons are squeezed towards 
the walls of our confinement. Different electron densities 
n^(r), n^(r) lead to different potentials, which, in turn, 
result in the difference in the conductance for two spin 
channels, Gj, (Fig.|3| left panel). 

The degree of spin polarization in open dots was 
probed by Folk et al. |7j who concluded that the dots 
were spin degenerate at low field. At the same time, 
the results of Rcf. H suggested that the spin degener- 
acy in open dots is lifted. Our findings reported above 
seem to reconcile these apparently different experimen- 
tal observations. Indeed, the dots studied in Ref. Q 
had relatively large electron density (n s ~2x 10 15 m~ 2 ), 
whereas the dots in Ref. had low electron density, 
(fi a < 7x 10 14 m~ 2 ), such that the observed features 
are consistent with our findings that spin polarization di- 
minishes as the electron density in the dot increases. 

During recent years a significant attention has been 
devoted to the study of the so-called "0.7-anomality" ob- 
served in quantum point contacts (QPCs) [(|. Though 
still a highly debated issue, this feature has been at- 
tributed by many researches to the effect of spin polar- 
ization that occurs at low densities when the exchange 
energy becomes more prominent compared to the kinetic 



4 




Conclusion. We have performed full quantum mechan- 
ical calculations of the magnetotransport in open quan- 
tum dots where the electron- and spin interactions have 
been incorporating within the local spin density approx- 
imation. At low electron densities the current is spin- 
polarized and electron density in the dot shows a strong 
spin polarization. As the electron density increases the 
spin polarization in the dot gradually diminishes. These 
findings are consistent with the existing experimental ob- 
servations. 
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FIG. 4: Average electron density, n], nf,, in the dot (left 
panel) and conductance (right panel) vs. magnetic field for 
spin-up and spin-down electrons calculated within LSDA for 
(a) low, (b) intermediate, (c) high density dots. B — 0T in 
(a), (b), and B = 0.01T in (c). T=2K. 



energy [ill fl3l lisj . Our present calculations show similar 
trends in the density dependent spin-polarization in open 
quantum dots, indicating the similar origin of the effect 
for these systems. Note that some recent studies ques- 
tion the existence of the spin-polarization in QPCs and 
attribute the "0.7-anomality" to different mechanisms, 
in particular, to the formation of one-dimension (ID) 
Wigncr crvstal[l5t. Such the interpretation relies on the 
notion of ID interacting system. In contrast, our findings 
strongly indicate that the spin-polarization at low elec- 
tron densities is a generic feature of both quasi-lD and 
quasi-OD systems (i.e. QPC and open dots). Because 
of this, we hope that our results will stimulate further 
experimental studies of the open quantum dots which 
might help to settle still unresolved issue concerning the 
spin polarization in the above systems. 

In our previous work Q , we investigated transport in 
similar dots using the mean-field Hubbard hamiltonian, 
H a = H n + J7rtf (r)nj(r), where on-site Hubbard con- 
stant U describes the Coulomb interaction between elec- 
trons of opposite spins. Within this model wc found that 
spin polarization in the dot takes place only when the 
Fermi energy hits the resonant states weakly coupled to 
the leads. In contrast, the present calculations based on 
the LSDA shows qualitatively different behavior, when 
the spin polarization is nearly constant in the broad in- 
terval of magnetic field. This difference strongly indi- 
cates that utilization of the simplified approaches that 
use phenomenological parameters and/or model Hamil- 
tonians might not be always reliable for theoretical pre- 
dictions as well as interpretations of the experiments. 
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